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Abstract 

Current-induced magnetic domain wall motion, induced by transfer of spin transfer effect due to 
exchange interaction, is expected to be useful for next generation high-density storages. We here 
show that efficient domain wall manipulation can be achieved by introduction of Rashba spin-orbit 
interaction, which induces spin precession of conduction electron and acts as an effective magnetic 
field. Its effect on domain wall motion depends on the wall configuration. We found that the effect 
is significant for Bloch wall with the hard axis along the current, since the effective field works 
as (3 or field- like term and removes the threshold current if in extrinsic pinning is absent. For 
Neel wall and Bloch wall with easy axis perpendicular to Rashba plane, the effective field induces 
a step motion of wall corresponding to a rotation of wall plane by the angle of approximately vr 
at current lower than intrinsic threshold. Rashba interaction would therefore be useful to assist 
efficient motion of domain walls at low current. 
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I. INTRODUCTION 



In recent years, magnetic-memory devices like hard disk drives have been utihzed for var- 
ious products such as portable music players and home appliances. With the diversification 
of the use of the devices, they deal with a much larger amount of information and need to be 
more miniaturized and have higher capacity. However, it is believed that the conventional 
magnetic-memory devices will reach the limit of downsizing and high capacity in the near 
future. The reason is because the magnetization so far is controlled by a magnetic field 
and the more miniaturized and the higher capacity they get, the higher necessary field they 
will need. Therefore, a new method which replaces the magnetic field is required. One of 
the methods expected is the control by electric current, i.e. current-induced magnetization 
reversal. The current- induced magnetization reversal was pointed out by Berger 
and Slonczeski. ^ The torque is caused by the s-d exchange interaction via spin transfer 
effect arising from conservation of spin angular momentum between the current and the 
magnetization. 

To realize devices using the current- induced magnetization reversal, it is absolutely es- 
sential to reduce the current necessary. In this paper, we demonstrate that the problem can 
be solved by introduction of Rashba spin-orbit interaction, A spin precession induced by 
Rashba interaction is expected to affect the spin transfer torque mechanism, and to lead to 
efficient control of magnetization. The interaction has been proposed on two-dimensional 
electron systems realized at the interface of semiconductors, but is now known to arise quite 

n 

generally when inversion symmetry is broken. [6| For instance, significant Rashba effect can 
arise on surface of heavy metals. [7|, l8| Our results therefore would apply to systems of mag- 
netic semiconductors 9| with gate for Rashba interaction attached, or on magnetic metallic 
thin films with heavy ions doped. 

Recent theoretical studies revealed that current- driven domain wall motion is significantly 
affected by spin relaxation process. Actually, spin relaxation triggers a torque perpendicular 



to the adiabatic spin transfer torque and this torque, /3 torque acts as an effective 



force on domain wall, which deletes intrinsic pinning effects ll|, ll2|. Microscopic analysis of 
spin relaxation was done in the case of spin flip scattering by random impurity spins and 
value of /3 was found to be of similar magnitude with (but not necessarily equal to) the spin 



relaxation contribution to Gilbert damping parameter a 



14 



16| . On the other hand. 
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role of spin-orbit interaction, which also causes spin relaxation, on the wall dynamics has not 
been much studied so far. Spin-orbit interaction in magnetic semiconductors was recently 
studied based on Kohn-Luttinger Hamiltonian, and large enhancement of wall velocity due 
to spin-orbit interaction was found [l^j. The Rashba spin-orbit interaction we are going to 
study here turns out also to assist wall motion to a large extent. In this paper, we study 
current-driven domain wall motion in the presence of Rashba interaction. Rashba interaction 
acts on x-y plane and electric current is applied along x axis. The spin system we consider is 

lich is justified if easy axis 



IS 



201 ]. We consider three 



with easy and hard axes. The wall is treated as planar and rigid, w : 
energy gain, is larger than hard axis anisotropy energy, K\ jl8. 
cases where easy axis direction (we call 77) is x, y and z. The wall structure in these three 
cases are Neel for rj = x and Bloch wall for t] = y,z. We will see that dynamics of Neel wall 
and Bloch wall with rj = z (called Bloch(z)) are essentially the same, while Bloch with t] = y 
(Bloch(y)) is different. The effect of current is calculated for each anisotropy configuration 
using gauge transformation in spin space, assuming adiabatic limit. The Rashba interaction 
is treated perturbatively to the second order. Expansion with respect to Rashba interaction 
is justified here since electron has large spin polarization A. (This is in contrast to the spin 
Hall case without polarization, where perturbative treatment is not allowed j2l|.) Below, 
we derive the effective Hamiltonian for local spin, derive the equation of motion for domain 
wall, and then discuss wall dynamics solving the equation. 

II. MODEL AND METHOD 
A. Local spin 

Domain wall is described by a Lagrangian of local spins given by 

-^hS {cos 9 -1)-H,, (1) 

where local spin direction is expressed by polar coordinates {6,4>)- The first term describes 
time evolution of local spin (spin Berry phase term). The Hamiltonian of local spin we 
consider is given as follows; 
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The first term is ferromagnetic excfiange interaction between local spins, the second and 
third terms are magnetic anisotropy energies, is the easy axis component of spin and 
K{> 0) is the corresponding anisotropic energy, and S_[_ and i^^ are hard axis ones. 

In the absence of Rashba interaction, we can choose space coordinate and spin coordinate 
independently. When Rashba spin-orbit interaction is switched, spin and space coordinates 
correlate each other, and the spin torque and gauge field depend on the choice of magnetic 
easy axis, which we call r] axis. We consider three different cases with easy axis t] = x,y 
and z. Domain wall structure then becomes Neel and Bloch walls. We will call domain wall 
configuration with the easy axis label, such as Bloch(z) for a Bloch wall with rj = z (see 
Table. ^. 





Easy (??) axis (K) 


Hard axis {K±) 


Bloch(z) 


z 


y 


Neel(x) 


X 


z 


Bloch (y) 


y 


X 



TABLE I: Configuration of magnetic anisotropy and corresponding domain wall structure. 

In this paper Rashba interaction acts within xy-plane and current is applied always along 
X. Instead, our definition of polar angle in spin space depends on the wall configuration. We 
define 6 as the angle measured from the easy axis (i.e., S\\ = cos 6'), and as the angle in 
the plane perpendicular to easy axis. The hard axis is given by </> = | and so the hard axis 
component is written as 5*^ = sin 6* sin 0. Note therefore that (6*, 0) for Neel and Bloch(y) 
walls below are different from standard definition measured from z-a.xis. Lagrangian Lg is 
thus written in terms of polar angle as 

hS(f) (cos - 1) - - {{Vef + sin^ ^( V0)2) - ^ sin2 6 - sin^ 6 sin^ , 

(3) 

for any type of walls considered here. 
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B. Conduction electron 



The Hamiltonian of the conduction electron is given by the following four terms. The 
first term is the free electron part, represented as 



X] 



c(a;), 



(4) 



where c and are electron annihilation and creation operators, e-p is the Fermi energy and 
m is the effective mass. The second is the exchange interaction between electron and local 
spin, 



ifcx = j d^xS ■ (c\x)cTc{x)) , 



(5) 



where A is the magnitude of the interaction, S is the local spin vector and cr are Pauli 
matrices. The third one, Rashba spin-orbit interaction, is given as 



Hso = - I d\c\x) (V.a^- V,a") c(x). 



(6) 



where V acts on both sides and A is the magnitude of Rashba interaction (the dimensions is 
[J ■ m]). When electric field is applied in the x-direction, the momentum, —ih (c^(x) Vxc(x)), 
grows, and Rashba interaction then changes the electron spin direction toward y-direction, 
inducing electron spin precession. 

Since we are interested in a response of local spins to applied current, we introduce the 
interaction with electric field, given by (neglecting 0{E'^)) 



EM 



-e J d^xApM ■ {c\x)vc{x)) 



(7) 



where Aem = ^c**^* is a U{1) gauge field, E is apphed electric field. Velocity operator is 



given as 



(8) 



The field E is spatially uniform but has finite frequency fi. This frequency is introduced for 
calculation purpose and is chosen as f2 = at the end of calculation, as is usually the case 
of linear response calculation. 
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C. Gauge transformation 



The exchange interaction, Hex, has in general off-diagonal components. Besides the local 
ground state of conduction electron varies at each lattice point if magnetization is non- 
uniform. In this case, local gauge transformation in spin space 3, 23] which diagonalizes 
the exchange interaction is useful. The spatial change of local spin is then represented by a 
gauge field, which is proportional to spatial spin variation, d^S. We consider the case when 
local spin profile is slowly varying (called the adiabatic limit), and then gauge field is small. 
We thus examine only the first-order contribution of gauge field. A new electron operator 



,a, after gauge transformation is defined by use of 2 x 2 unitary matrix U as 

c(x) = U{x)a{x), 



(9) 



where electron operators here have two spin components like c = (c+, c_). The matrix U is 
expressed using Pauli matrix as 



U{x) 



mix) 



(10) 



where m(x) is a vector which characterize the gauge transformation. We denote the 
spin easy axis as rj (e.g., t] = x for Neel wall). The gauge transformation is defined so 
that conduction electron spins are polarized along magnetic easy axis rj, i.e. to satisfy 
■ (c^(x)crc(x)) = a^(x)(T^a(x). (The transformation thus differs for different wall config- 
uration.) The Hamiltonian given by Eqs. (|1])([5])(I6]) is written in a-electron representation 
as 



He + He 



H^, 



d^x 





2m 



m 



where gauge field is given as 



A., 



Aa" - zA(aW^ - a^Vy) 
+ 2zA [mym'^a^a^Vc.a - m^'m^a^a'^Vya] 



-iU^V^U = (m X d^m) ■ a- = A^a'^ . 



The electromagnetic interaction after gauge transformation is given as 



H 



EM 



d'^x- 
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-ihVx ^ , A 

+ -Ax - -ay 

m m n 



(11) 



(12) 



(13) 
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where o"^ = Wcr^U = 2m^ (m ■ cr) — and the electric field is applied in x-direction. 
The total electron Hamiltonian is therefore obtained as 



H = (Tx 



2m Vt \ m ih 



m 



in ' 



—Ay - (m ■ a) 

m h 



(14) 



D. Effective Hamiltonian 



The current-induced part of the effective Hamiltonian for local spin is directly obtained 
from Eq. f|T^ as 



eff 



(fx 



hAl3l^ + \n{Al-Al) 



2m 



A [mymf^jl^ - m'^m'^jly] 



iVt 



m n 
(15) 



m^m s 



where electron properties are represented by the following expectation values, 

n{x) = (a^(x)a(x)) , 



-^(at(x)V,a-a(x: 



(16) 
(17) 

(18) 



Here ( ) denotes expectation value evaluated using non-perturbed Hamiltonian Ha defined 

as 



d xa'(x) 



V - Aa" - iXia^W^ - a^'VJ 

2m 



a(x) 



and including to linear order the effect of current, 



-a' 



Current-induced part of the effective Lagrangian is given by 



a{x). 



(19) 



(20) 



(21) 



In calculating the expectation values, Rashba spin-orbit interaction is treated perturba- 



tively to the second-order, A^. This approximation correspond to assuming ^ 1 with 



7 



Fermi wave vector kp. This expansion with respect to A is justified by the presence of A, in 



contrast to the non-perturbative nature of unpolarized Rashba system[2]|. We also assume 
that the effect of impurities is weak and the electron lifetime is long, i.e. ^ <^ 1. In this 
case, the term in Eq. ( IT^ including electric field turns out to be small by a factor of ^ 
compared with dominant contributions. We will thus evaluate the effective Hamiltonian 
given by 

2m 



eff 



-A K^^^t 



m m 



■'s,y\ 



(22) 



III. BLOCH(Z) CASE 

In this section, we derive the effective Hamiltonian for the anisotropy configuration of 
Bloch(z) type, namely, magnetic easy axis is in the 2;-direction and hard axis is in the y- 
direction. In this case, domain wall configuration is Bloch(z) wall shown in Fig. [TJ Polar 



FIG. 1: Bloch(z) wall configuration with 
easy axis along z direction and hard axis 
along y axis. The spin of incoming electron 
is schematically shown. 




FIG. 2: Polar coordinates of local spin in the 
Bloch(z) case. 



coordinates are defined by (as in Fig. [2l) 



S{x, t) = S{sm 9 cos 0, sin 9 sin 0, cos ( 



(23) 



In the case of Bloch(z)-type anisotropy, the gauge transformation is defined by 
S{W(tU) = Sa^ . This is satisfied by choosing 



m = sm - sm 0, cos - , sm - cos i 
^2 2 2 



(24) 
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Gauge field is then given as 







1 

~ 2 













^ (?^6'sin0 — sin6'cos</)(9^(/'^ 



d^O cos ( 



sin 9 sin < 



V 



(l - cos9)d^(, 



(25) 



The unperturbed Hamiltonian, given by choosing r] = z in Eq. f|T9|) . reads 

^2 



X 



2m 



a{x) 



(26) 



This Hamiltonian i/a has the off-diagonal elements in the spin space due to Rashba inter- 
action, and so we will diagonalize it using a unitary transformation in momentum space, 
dlrl, where dl is a new creation operator and is 2 x 2 unitary matrix. The 



Hamiltonian Ha in momentum space reads 



iXk^ Ek — A 



Tkdk 



(27) 



where k± = kx i: iky. Diagonalization of Ha is done by chosen the unitary matrix Tk as. 



Ak iXk. 



-Au. 



(28) 



where A^ = A + Z^., = vA^ + A^A?^. The result is 




(29) 



After diagonalization, the energy of conduction electrons is spin split as 

hk^ 



2m 



Ep — crZ{k), 



(30) 



where a = (+, — ) is spin polarization. 

Now, let us estimate expectation values. The electron density n (Fig. [3]) is given as 

1 < 



n 



-|- lim > — Tr 



ci^o ^-^ fl m 
+ lim V ^ATr 



1 < 



(31) 
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FIG. 3: Diagramatic representation of electron density (n), spin density (s") and spin current (jf J 
at the linear order in the applied electric field, E. Solid line represents electron Green function 
with Rashba interaction included. Vertex denoted by x is 1, fi" and /cjO"", for n, s°' and 
respectively. 



where tr is a trace in spin space and 5"° = Tla°'Tk. Keldysh Green function 24, 25| with 



Rashba interaction included, G^^, is defined as 




(32) 



where component g^f^ given as 



(33) 



/(c<j) = -^^u^ being Fermi distribution function (/3 is inverse temperature). Summation 
over k is carried out in two- dimensions by replacing by energy integration, Ylk ~ I 
Products of Keldysh Green functions is calculated using relations 



(34) 



where ^ {G^. ^) is retarded (advanced) Green function, and 



V^k,uj^k',u)+Q.\ —^k,u^k',u)+n + ^ k,u)^ k' 



— ^k,l^^k',i 

=Gk,ujGk',u;+n (/(^ + ^) - f{uj)) - Gl^^Gl,^^_^^ + Gl^^Gl,^^_^^. 



(35) 



Expanding with respect to A to the second order, we obtain the density as 



m 

n = -^€p + 



We see that there is no effect from the applied electric field here. 



(36) 
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The electron spin density s'' is similarly calculated as 
s7 = _^,Tr b^Gfc, 



+ lim > — Tr 



\ — ^ (^Ex A 
+ lim > — — — ir 



(37) 



The result is 



meE^rX 

— 

Txh^ 

niA 



We see that the electric field induces perpendicular components and s^, but 
since ^ 1- We will thus approximate |s^| ~ 0. The 2;-component of spin in Eq. 
the adiabatic contribution, which is not affected by applied field. 
Spin current, j"^, is estimated using 

< 



(38) 

(EHD is 



k,LO 



— lim > k., — Tr 



iVtm 



'^°'Gk,ioGk,Lo+n 



lim k,, ^ j!^ Tr 



0'°'Gk,ujO'^Gk,uj+Q 



(39) 



as 



■Js,y 


m 




m 


J s,x 






eExT 


J s,x 





A 



(40) 



We see that spin current are generated in x and y direction by the Rashba interaction 
without electric field and that j^ y is equal to —jj^^. This is due to the symmetry of Rashba 
spin orbit interaction. In contrast, ^ is induced by applied electric field. We define the 
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current density (divided by e) and electron density without spin-orbit interaction as {a is 
Boltzmann conductivity) 



J 



eExTS-p a 



Txh 



(41) 



and 



no 



In term of these parameters, the above result reads 



(42) 



n =nQ 



1 + 



'^^-"-•)HO.-^i,«o|) 



(jy f ^ 



X 
h 



(43) 



Other components of spin current vanish. 

The effective Hamiltonian for Bloch(z) case is therefore obtained from Eq. fl22|) as 



cff 



(fx 
Am 



h 



l-cos^) {d^<t>)jl 



[jIx sin ^ sin </) + (1 - cos 9) [jI^ - sin (pjly)] 



(44) 



We see that applied current induces j^^ (Eq. fj^3|l ) and this induces when coupled with 
Rashba interaction an effective magnetic field in ^/-direction as indicated by the second 
term. The first term of Eq. (jHj) represents standard spin transfer torque (with current 
modified by Rashba interaction). The third term is independent of applied current and is a 
modification of magnetic anisotropy by Rashba interaction. 



IV. NEEL(X) CASE 

In this section, we consider a case of Neel wall realized when magnetic easy axis and 
hard axis are in the x and z direction, respectively (Fig. H]). Polar coordinates are defined 
differently from standard definition as (Fig. [5]) 

S = S{cos 6, sin 6 cos sin 6 sin (p) . (45) 
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FIG. 4: Neel(x) wall configuration. 



FIG. 5: Definition of polar coordinates for 
Neel(x) case. 



Derivation of effective Hamiltonian is done similarly to Bloch(z) case. Difference is in 
definition of gauge transformation, SiJJ'^crU) = Sa^. Vector m is accordingly chosen as 



e e 



e 



m = [cos -, sm - cos (p, sm - sm i 



(46) 



and is given as 



/ 



A,. 



(1 — cos 6)d^(j) 



— sin (f)d^9 — sin 9 cos 09^^ 
y d^9 cos — sin ^ sin 09^0 j 



(47) 



Hamiltonian Ha is also different from Bloch(z) case since uniform spin polarization is now 
along x-direction. It is given as 



where b-i 



-b+ Sk 

(A =1= i\k±). The diagonalization of Ha is carried out as aj[. = d\Tl where 

1 6_ 



(48) 



T = 



where 



A2 + 2\/\ky + X^k"^. 



(49) 



(50) 



The Hamiltonian after diagonalization reads 



' eu + Zu 



(51) 
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The expectation values are calculated similarly to Bloch(z) case and we obtain 



n =nQ (52) 




where other components of spin and spin current vanish. 

The current-induced effective Hamiltonian is then obtained as 

HeS= J rf^(l-COS^) 

- -A [jI, sin cos (/) + (1 - cos 9) {jI, cos" (j) - jI^) ] . (55) 

The applied current coupled with spin-orbit interaction induces an effective magnetic field 
along hard, i.e., y axis (second term of Eq. fl55i) ). as in the Bloch(z) case. 



V. BLOCH(Y) CASE 




FIG. 7: Polar coordinates for Bloch(y) con- 

FIG. 6: Bloch(y) wall configuration. 

figuration. 

Finally, we consider a case of magnetic easy axis and hard axis along y and x direction, 
respectively. Domain wall in this case is Bloch(y) wall (Fig. [6]). Polar coordinates are define 
as (Fig. 

S = S{sm 9 sin 0, cos 6^, sin 6^ cos 0) . (56) 

Gauge transformation is given by S{WaU) = Scr^, and vector m consequently becomes 

9 9 9 / 

m = (sin - sin 0, cos -, sin - cos 0), (57) 
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and gauge field is 



Hamiltonian Ha is given as 



^ 3^9 cos — sin sin 09^0 ^ 



'1 — cos 6)8^ 



— sin (f)d^9 — sin 6' cos 0(9^ 



(58) 



-ib+ Ek 



(59) 



where 6± = (A =i= Xk±) (is different from the Neel case). Hamiltonian Ha is diagonalized as 



£k — Zk 



by defining a\, = d^Tl with 



-Zk ib- 



V2Zk \-ib+ Zk 



where 



(60) 



(61) 



Zk = - 2XAk^ + X^k^. 



(62) 



Expectation values which appears in the effective Hamiltonian are estimated as 

3m 



n =nQ 



X' 



A mA^ 



X 



XJ 
A mA^ 



2mA2 



Eph? Ah? 



X mX^ 



h At? 



J 



The effective Hamiltonian is therefore obtained as 

'1 



H 



cfT 



mX 



l-cos^) (9^0) [hil^ + Xn] 
[(l + cose)jJ,,-(l-cose)j-,sin2 



(63) 



(64) 
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VI. ANALYSIS OF EFFECTIVE HAMILTONIAN 



The effective Hamiltonian representing the effect of current obtained above is summarized 
as follows. 
Bloch(z) case: 



H, 



cS 



+ 



2 ^ 
no (1 



cos 6*) ((9^0) 
Cz) (1 + sin (J 



A 

Bp 



m 



2Ah? 



J 



AmA 



(65) 



Neel(x) case: 



eff 



(fx 



^(1 

2 ^ 



cos I 



A AmA 

[d^(p) — J —CyJ 



.) no (l + cos^ (f) 



(66) 



Bloch(y) case: 



(fx 



h 



cos 9) (9^0) 



LA m 

5f £^F^^ 



A^ + 



2m 
An2' 



J 



AmA 



a2 r 
+-^A no [ey + 



sm 



(67) 



Here Cq is local spin component in a-direction. The first term of each effective Hamiltonian 
represents spin transfer torque, which is enhanced by Rashba interaction in Bloch(z) and 
Bloch(y) cases at the second order in spin-orbit interaction. The enhancement is thus small 
but independent of sign of A. The second term of the effective Hamiltonian, CyX-^J, indicates 
that the effective magnetic field arises in the ^/-direction in all three cases. This is a result 
of spin current induced by Rashba interaction and applied current. The last term of each 
Hamiltonian exists without electric field and thus represents magnetic anisotropy modified 
by Rashba interaction. This change of magnetic anisotropy is at the second-order in A and 
is small. Furthermore, in reality, this static contribution should be contained already in the 
anisotropy parameters K, K±. We therefore do not take it into account in the following 
analysis. 

The effects of Rashba interaction are summarized in Table. [TTl 
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Bloch(z) 


Neel(x) Bloch(y) 


Enhancement of spin transfer torque 




X 0(A2) 


Effective magnetic field 


0(A) 


0(A) 0(A) 


Change of magnetic anisotropy 


0(A2) 


0(A2) 0(A2) 



TABLE II: Summary of effects of Rashba interaction for three configurations represented by the 
order of the effects in A. Symbol x denotes the absence of the effect. 

VII. EQUATION OF MOTION OF DOMAIN WALL 



Let us discuss how Rashba spin-orbit interaction affects current-induced domain wall 
motion based on the effective Hamiltonian we derived. For this Lagrangian formalization 
is convenient. We consider rigid planar domain wall. The Lagrangian for domain wall is 
then obtained as L = Xcf) — {Hs + Hcs)\x,(f> where X is domain wall position, (p is local spin 
angle in magnetic easy plane and {Hs + -f^cfr)U,(^ is effective spin Hamiltonian evaluated for 
domain wall configuration [isl . Explicitly, domain wall Lagrangian is given as follows. 
Bloch(z) wall: 



L 



(f)X — sin^ I 
- s'X^ - 4 



ITS 



y 



+ cos (pep - 
A 



2 



X cos( 



sm( 



A ~ 

J -|- vrA sin0 — J. 



(68) 



Neel(x) wall: 



(pX - sin^ 



s''X(j) + 



Bloch(y) wall: 



(pX - sin^ 



A 



X sin 6 + cos i 



.X r 



!)— J + 7rAcos0— J. (69) 



J - 2A 



~X{t) A ~ 



Here we introduced following dimensionless parameters , 

X r X 



J. 



L 
J 



hNSvc 
J 



A 



X = 

mi 



X 



A. 
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where is a number of local spin in the domain wall and ^ is domain wall thickness given 
as £ = \ j^- (Dimensionless time is t = yt). Here a velocity Vc = ^on^ correspond to drift 



velocity of electron at intrinsic threshold current without Rashba interaction [l8]. 



Equations of motion for domain wall is obtained taking account of dissipation a.s,{Q 



d dL dL 
dtd^ ~ dQ 



' dQ 



(70) 



where Ws is dimensionless dissipation function written as Ws 

In deriving the equation of motion, we neglect contribution of conduction electron density 
{s'^ <^ 1) since they turns out to be small in actual situations. The equation of motion is 
obtained as follows. 
Bloch(z) wall: 



X — ac 



- sin 20 



ef 2AmP 



A , 

J — irX cos</) — J, 



(71) 



+ aX = 0, 



(72) 



where time evolution of (p reduces to a single equation of 

■A 



a 



Neel(x) wall: 



1 + a2 



sin 20 



Ef 2AmP 



A ~ 

J — nX COS0 — J 
Ef 



A ~ ~ A ~ 
X — acj) = sin 20 H J + ttA sin0 — J, 

Ef £f 



(73) 



(74) 



+ aX =0, 



(75) 



which result in 



Bloch(y) wall: 



a 



A ~ ~ A ~ 
sin 20 H J + TxX sin — J 

Ef £f 



X — ad) = sin '. 



A 



h^x^ 



2h^X^ 



Ef tjiPef mPA 



J, 



(76) 



(77) 
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2AA ~ 
+ aX = J, 



(78) 



and the time evolution of ch is obtained as 



a 



sm . 



A 



J 



2h\^ X 2~XA ~ 
-no sm 20— H J 



m£vr 



a Ef 



(79) 



We here see a large difference between Bloch(y) and other cases. In fact, the angle in 
Bloch(y) wall is directly driven by current and Rashba interaction as indicated by the right- 
side of Eq. (178|) . Such effect of current has been known as f3 terms in the case of electron 



spin relaxation due to random spin 
then given by 



10 



26| . In the present Rashba case, the parameter j3 is 



2AA 



(80) 



for Bloch(y) (/5 = for Bloch(z) and Neel(x) walls). This result can be explained by 
noting that (3 terms is effectively equivalent to an external magnetic magnetic field and 
that Rashba interaction induces an effective magnetic field, which coinsides with the easy 
axis for Bloch(y) wall. (Eq. (1671) . We will see below that this coefficient j3 is quite large 
even assuming standard semiconducting systems and reduces much the threshold current 
and enhances the wall velocity. The effect would be even stronger if the systems has giant 
Rashba effect as realized in metallic surfaces 81] . 

The effect of effective magnetic field due to Rashba interaction in Bloch(z) and Neel(x) 
cases is to induce anisotropy within 0-plane as seen as cos0 and sin0 terms in Eqs. (ITTll 
and (171|) . respectively. This anisotropy energy turns out to drive stepwise wall motion at 
low current. 

We also see from Eqs. (17T|) and (1771) that spin transfer torque effect (represented by the 
second term of right-hand side) is enhanced by Rashba interaction, at the second oder of 
A in Bloch(z) and Bloch(y) cases, but not in Neel(x) case. Numerically, these second order 
effects are negligibly small as we will show below. 
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VIII. DOMAIN WALL DYNAMICS 



We first note that domain wall velocity X is closely related to <p 18|. In fact, from Eqs. 
(IZ2]),([75]),([7S]), we see that 



and 



X = --(p (Bloch(z) and Neel(x) wall), 
a 



X = --<j) - ^ J (Bloch(y) wall). 



(81) 



(82) 



Here, we show numerical results based on the equation of motion. Parameter are chosen 



as 



a = 5.65 X 10 ^°m, m = 0.05me 
Vc = 163m/s, ep = SOmeV, A = 40meV, 



me = 9.11 X 10~^^kg, 
' = 50nm, a = 0.01. 



to simulate actual semiconductor systems [9|. (We thus have 



^2 



0.007, and so the 



second order contribution from Rashba interaction is very small like A 



2 ^2 



0.7 X 10"^ 



for A = 0.1.) We calculated the domain wall position X{t) after current is applied at t = 0. 
Fig. [8] shows the wall position at t = l/im after current is applied. (Current is normalized by 
Vc, the electron drift velocity at intrinsic threshold.) The case without Rashba interaction is 
shown as solid line. We see here the intrinsic pinning due to hard axis anisotropy 18|] , since 
we do not consider jS term of non-Rashba origin. Wall motion in the presence of Rashba 
interaction with A = 0.1 is plotted by lines marked by x for Bloch(z) and Neel(x) and * for 
Bloch(y) walls. Bloch(z) and Neel(x) walls behaves essentially the same. 



A. Bloch(y) wall 

We immediately see that Rashba interaction affects Bloch(y) wall drastically, resulting 
in vanishing of threshold current and very high velocity. This is due to a large effective 
f] term induced by Rashba interaction, Eq. fl80l) . (Note that we do not consider extrinsic 
pinning.) For the present parameters, its ratio to a is given as £ = ^, and so very large 
value of ~ 10 can be realized in actual experiment with A ~ 0.1 and a = 0.01. Terminal 
wall velocity of Bloch(y) wall is plotted for different values of A in Fig. O The behavior is 
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FIG. 8: Wall position at t = Ifis as function of dimensionless current, J = J/v^. Solid line 
represents the case A = 0, and dashed and dotted lines represent the case A = 0.1 for Bloch(z), 
Neel(x) and Bloch(y) walls, respectively. 



consistent with previous studies including (3 term ll|, ll2l] indicating that the present motion 
of Bloch(y) is governed by (3 term induced by Rashba interaction. It is natural since (3 term 
is linear in A Eq. f lHOj) while other effects of Rashba interaction enter at the second order (see 
Eq. (177j) ). Within our analysis, which neglects the extrinsic pinning, threshold current of 
Bloch(y) is zero due to (3 term. In reality, however, finite threshold will appear from pinning 
potential jlll, [l^. Due to the (3 term, we have obtained very large wall velocity (like lOOm/s) 
at very small current (like 10% of intrinsic threshold current). In reality, Gilbert damping 



a) is also modified by Rashba interaction, which might slow the ve 



V oc /9/a|12l. The modification of a is, however, second order in A 
hence the high wall velocity would remain unchanged. 
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o city V according to 
28| and thus small. 
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FIG. 9: Plot of wall velocity as function is current for A = 0,0.001,0.005,0.01 and 0.1. Damping 
constant is a = 0.01. 

B. Bloch(z) and Neel(x) walls 

Behaviors of Bloch(z) and Neel(x) walls are essentially the same for present values of 
parameters, and are distinct from Bloch(y) wall. These walls have a plateau in v-J curve 
for 0.8 < J < 1.3 as seen in Fig. [HI This plateau is due to a step motion of wall induced by 
the anisotropy field within 0-plane arising from Rashba interaction. This stepwise motion 
is induced above threshold current of ~ 0.8, and the distance the wall moves is about 22^m 
regardless of current density (for 0.8 < J < 1.3). 

Let us see the mechanism of the plateau in detail. Dynamics near intrinsic pinning 
threshold is described by a potential for obtained from the Lagrangian (Eq. (l68l) and 

Eq. dSni)) as follows: 
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Bloch(z): 



y(0)=sin20+^ 



Ef 2AmP 



TT ~ A ~ 

J - —Asm 0— J. 



(83) 



Neel(x): 



A ~ TT r 



A = 



=sin 0+777; — J — Trp;X cos (j) — J 



(84) 



25 25 ef 

These potentials are plotted in Fig. [TOl In the absence of current, the potential is solely 



(1) 







2x wr 




FIG. 10: Potential for of domain wall. Figure (1) is the potential without current and Rashba 
interaction, namely potential due to the hard axis anisotropy energy. Energy minimums are at 
(j) = 0, TT, • • • . Figure (2) includes the effect of applied current below threshold, which tilts the 
potential (J < here). (The intrinsic threshold at J = 1 corresponds to the current where local 
minimum near ~ disappears.) The last figure (3) shows the potential in the presence of Rashba 
interaction, which lowers (lifts) the energy barrier around ~ § for Neel(x) shown by dashed line 
(Bloch(z) shown by solid line). 

by magnetic anisotropy, and energy minimum is at = (Fig. [TOT l)). When current is 
applied, the potential tilts (Fig. [10] (2)). When local minimum near = vanishes, the 
wall starts to move, and this gives the intrinsic threshold current of Jc = 1 in the absence 
of Rashba interaction, [l^ When Rashba interaction is switched, vr-periodicity of potential 
is broken due to the last terms in Eqs. fl83l) (!8^. The way of deformation depends on 
Bloch(z) and Neel(x) cases. Let us consider a Neel(x) wall case where the local energy 
barrier around vr ~ ^ is lowered by Rashba interaction. The variable (and wall) starts 
to move at threshold current (where the local minimum disappears), which is lowered by 
Rashba interaction. But this motion stops if current is not large enough since is trapped 
by the next local minimum near ~ vr, which has a large energy barrier around tt ~ 
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Thus in this regime, hops roughly by the amount A0 = 0(vr) and this corresponds by 
Eq. fjHTl) to a wall shift of AX = ^A0. This distance corresponds to the distance ~ 20/i m 
seen in Fig. [8l If current is sufficiently large to remove the second local minimum, the wall 
motion continues and terminal velocity grows (for J > 1.3 in Fig. [S])- This is the mechanism 
of plateau. 

An interesting consequence of the step motion is an asymmetric ratchet motion. We 
consider a Neel wall initially at </> = (we choose A > 0). When we apply a negative 
current, J < 0, the local minimum at = is lifted and the energy barrier for right 
direction is lowered as we saw in Fig. [101 In contrast, when current is positive, J > 0, the 
local minimum is lowered (oc —\ J by Eq. (18^ ) and then the effective energy barrier in 
the left direction becomes higher. Therefore, the threshold current for step motion differs 
by amount A J ~ A for the two current directions (starting from fix 0), and thus the wall 
behaves as a ratchet moving only in one way if current is small enough. These features are 
common for Bloch(z) wall case if = is replaced by = vr. 

The plateau region has finite initial velocity (but zero terminal velocity). Figure fTTl shows 
the initial and terminal velocities for Bloch(z) and Neel(x) walls (the velocity is the same for 
two walls). We see that even plateau region shows a high velocity comparable to terminal 
velocity above J = 1. For device application, motion over a distance of 20/im is large enough 
and so the plateau region would be quite useful for low current switching. 

IX. CONCLUSION 

We have theoretically calculated the effect of the Rashba spin-orbit interaction on the spin 
transfer torque. The effective Hamiltonian of local spin under current was calculated using 
gauge transformation, and the equation of motion for domain wall was derived. We consid- 
ered three cases with different magnetic easy and hard axes, where domain wall structures 
realized are called Bloch(z), Neel(x) and Bloch(y) walls. We found there are three influences 
of Rashba spin-orbit interaction, namely, inducing effective magnetic fleld, increasing spin 
transfer torque and modiflcation of magnetic anisotropy. The major effect is that of effective 
magnetic fleld, which arises at the flrst-order in Rashba interaction. Applying voltage in 
x-direction, the effective magnetic fleld is induced in y-direction via Rashba interaction. In 
case of Bloch(y) wall, where y-direction is the easy axis direction, we showed this fleld acts 
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FIG. 11: The initial and terminal velocity as function of dimensionless current in the case of 
Bloch(z) and Neel(x) wall. The velocity is the same for the two walls. 

as a force which pushes the wall, or in other words, effective (3 term arises. Threshold current 
thus vanishes. The value of (3 is large even if evaluated for common semiconductor systems, 
and wall velocity is enhanced greatly. In contrast, in the cases of Bloch(z) and Neel(x), 
the effective field is perpendicular the easy axis {z and x directions, respectively), and step 
motion of wall over a distance of AX ~ O(^) arises at low current regime, corresponding 
to a change of the angle out of easy plane, A0 ~ 0{7r). The current necessary for this 
step motion is lower than the case without Rashba interaction (by 20% at A = 0.1). The 
initial velocity of step motion is high enough (the same order as steady motion slightly above 
intrinsic threshold). Wall motion in the step motion regime is asymmetric with respect to 
current direction, i.e., wall behaves as a ratchet. 

Other effects by Rashba interaction, modification of spin transfer torque and anisotropy, 
appears at second order, A^, and are negligibly small. Change of spin transfer is due to the 
change of effective electron spin polarization by Rashba interaction. 

Rashba interaction arises quite generally when inversion symmetry is broken, e.g., on 
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surface of metals doped with heavy ions[7|, l8|]. Such systems would be suitable to realize 
quite high wall velocity at very small current we have predicted for Bloch(y). 
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